Interaction-Induced Quantum Dephasing in Mesoscopic Rings 
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Combining nonperturbative techniques with Monte Carlo simulations we demonstrate that quan- 
tum coherence effects for a particle on a ring are suppressed beyond a finite length L v even at zero 
temperature if the particle is coupled to a diffusive electron gas by means of long range Coulomb 
interaction. This length is consistent with L v derived from weak-localization-type of analysis. 
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Experiments Q strongly suggest that decoherence 
time Tip of electrons in disordered conductors saturates 
to a finite value in the low temperature limit. Theoret- 
ical arguments 0] attribute this saturation to the effect 
of electron-electron interactions. This work has initiated 
intensive theoretical debates, the present status of which 
is summarized in recent publications |s| . 

Although detailed calculations of t v for metallic con- 
ductors with disorder and interactions can be rather in- 
volved, the main cause for electron dephasing is quite 
clear without going into unnecessary complications: It 
is the electron interaction with the fluctuating quan- 
tum electromagnetic field produced by other electrons. 
A question of fundamental importance here is to under- 
stand whether such interaction can dephase at T — > 0, 
i.e. if the whole interacting system is in its true quantum 
mechanical ground state. 

It is obvious that for a system in thermodynamic equi- 
librium quantum dephasing cannot be associated with 
any kind of decay or relaxation in real time, rather it 
manifests itself via suppression of off-diagonal elements 
of the electron density matrix beyond a certain length 
Lip. Provided there exists nonzero electron dephasing 
due to its interaction with quantum environment even 
at T — > 0, this dephasing length should stay finite 
down to zero temperature. This, in turn, would imply 
that effects sensitive to quantum coherence, such as, e.g., 
persistent currents (PC) and Aharonov-Bohm (AB) os- 
cillations in mesoscopic rings, should be suppressed by 
interactions at any temperature including T = if the 
perimeter of the ring exceeds L v . 

Recently the authors M demonstrated that the am- 
plitude of PC can be reduced by interactions even ex- 
actly in the ground state, and this result was interpreted 
as a signature of suppression of quantum coherence at 
T = 0. Other authors || argued that this reduction of 
PC is merely a renormalization effect which has nothing 
to do with dephasing. Very recently Guinea || found 
that AB oscillations for a quantum particle on a ring 
interacting with Caldeira-Leggett bath of oscillators are 
exponentially reduced beyond some length L v which is 



set by interactions and remains finite down to zero tem- 
perature. Similar result was also obtained earlier from 
a real-time quasiclassical analysis J?J. It should also be 
noted that the model HQ is exactly mapped onto that 
of the so-called single electron box, the effective charging 
energy of which is well known to reduce exponentially at 
large conductances M O. This reduction immediately 
translates into a finite dephasing length L v at zero tem- 
perature found in On the other hand, no finite L v 
at T — was found in some other models |^,|l2| . 

The purpose of this Letter is to investigate quantum 
dephasing for a model with long range Coulomb interac- 
tions in the presence of disorder. Our main conclusion is 
that such interactions can strongly suppress PC even in 
thermodynamic equilibrium at T = 0. 

The model and effective action. We will consider a 
quantum particle with mass M on a ring with radius R 
threaded by external magnetic flux <& x . As in ref. Q it 
will be convenient to describe the particle position by a 
vector R = (Rcos9,Rsm9) and consider the angle 8 as 
a quantum variable. We also assume that the ring is em- 
bedded into an effective dissipative environment and the 
whole system is in equilibrium at a temperature T. 

Our first and standard step is to integrate out the en- 
vironmental degrees of freedom. After that the grand 
partition function of the system takes the form 
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The first term in the exponent takes care of the magnetic 
flux while the second term 
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is the particle action. Here we defined (3 = 1/T, Ec — 
1/(2MR 2 ) and <f> x = $ x /®o, where $ is the flux quan- 
tum. The term 5i nt describes the effect of interaction 
between the particle and the environment and, hence, 
depends both on the environment Hamiltonian and on 
the form of the interaction. We will assume that the 
particle has the electron charge e which interacts with 
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the fluctuating electromagnetic field V produced by the 
effective environment. For this model one finds 
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ln/exp i-i J dreV{r,6{T)) 
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Provided the fluctuating field V is sufficiently well de- 
scribed within the Gaussian approximation, we obtain 
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\t f dT'{V(T,e(T))V(r',0(r'))). (4) 
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The correlator (VV) can be expressed via the dielectric 
susceptibility of the environment e(u>,k) as: 

[VV> (2^k 2 e(i\u n \,k) 
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where w„ = 2?mT and X = R(t) — R(r'). In what 
follows we will model the environment by a 3d diffusive 
electron gas with l/e(w,fc) w (— iu> + Dk 2 )/ (iira), where 
a is the Drude conductivity of this gas and D = vfI/3 is 
the electron diffusion coefficient. Substituting l/e(i|o;| , k) 
into (@) and performing integrations we observe that 
the term ~ Dk 2 yields only an A-independent energy 
shift. Integrating the remaining contribution ~ \u} n \ over 
k < 1/Z in (||) we obtain the result oc min(l/X, 1/7) » 
l/VX 2 +l 2 . Defining X = 2iMn[(0(r) - 6»(r'))/2] and 
summing over cu n , from (Q) and (JsJ) we get 

o rfll [ P . . ^T^r)-^)) 

5i„ t = a / dr / cZt . 2 -— , (6) 

Jo Jo sin [t:T{t-t')\ 
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where a = i/(%k 2 F l 2 ) and r = R/l. The integral in eq. 
(^J) is understood as a principal value. The divergence 
at r = t' is then regularized in a standard manner by 
requiring K(0) = which explains the origin of the first 
term in (]?]). Below we will set 1/kp <§C / <C R imply- 
ing that interaction is effectively weak a <^ 1 while the 
parameter r = R/l is large r> 1. 

Perturbation theory. In the non-interacting limit 
S'int — * the partition function (0) is trivially evaluated, 
and for the flux-dependent part of the ground state en- 
ergy at —0.5 < 4> x < 0.5 one finds E (cj) x ) = Ecfil- 

Let us now take interaction into account. We first as- 
sume that interaction effects are weak, in which case it 
suffices to expand the partition function (|l|) to the first 
order in S'int ([!)■ Here it will be convenient for us to 
rewrite the function K in terms of the Fourier scries |q] 
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n(d(r) - 9(t')) 
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where a n ~ (2/wr) In ir/n) for 1 < n < r and a n « 
otherwise. After that the calculation reduces to a sim- 
ple Gaussian integration over the 9- variable in each term 



of the series. Performing this integration and summing 
over n one obtains Eq(4> x ) and finds a diamagnetic cur- 
rent I = (e/2n)(dE /d(f> x ) with the result 
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where —0.5 < a < 0.5. The last term in (g) repre- 
sents the first order correction to PC in the ring due 
to Coulomb interaction at T = 0. This correction is 
negative, i.e. for the problem in question interactions 
suppress PC even in the ground state. For <fi x <C 1 the 
interaction term in (|9|) is linear in <f> x and it is small in 
the parameter a Y] , ~ q <C 1 as compared to the 
non-interacting PC. A somewhat stronger effect of in- 
teractions ~ a In r was found in ref. Iq] within the renor- 
malization group (RG) analysis. Below we will show that 
already for small a <C 1 suppression of quantum coher- 
ence is actually much stronger than it could be expected 
both from eq. ([|) and the RG approach M . 

Nonperturbative effects. Let us first consider the 
regime of not very low temperatures. In this case the 
partition function (|l|) can be evaluated semiclassically. 
As in jlO 11| one can find the classical paths providing 
the minimum of the total action S = So + S'int for any 
winding number m. These are the so-called "straight 
line" paths 9 c \(t) = 2irmTT, describing rotations of the 
particle around the ring. Substituting the paths c \(t) 
into the sum (Q) one gets 
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In the limit T Ec / tt 2 it suffices to keep only the terms 
with m = 0,±1 and the flux-dependent part of the free 
energy F = — TlnZ can easily be evaluated. Then for 
PC I{T) = (e/2ir)(dF/d(j) x ) one obtains 



I = 2eTexp - 
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— Anar I sin(27n^ x ). 
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We observe that in addition to trivial reduction of PC 
due to thermal fluctuations there is also a T-independent 
term which yields exponential suppression of / provided 
the perimeter of the ring 2itR exceeds the value 



L v ~ l/a <~ l(kpiy 
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Eq. (|lj) defines the dephasing length for the problem in 
question. This length does not depend on temperature 
and is controlled by the effective interaction strength a. 

The above semiclassical analysis is justified at suffi- 
ciently high temperatures, whereas at lower T < Ec 



2 



fluctuations around the classical paths 9 c \(t) become im- 
portant and our treatment needs to be modified. Pro- 
vided the suppression of PC is sufficiently strong one can 
employ the instanton technique |P,[To|. 

Nontrivial saddle points 9(t) which describe quantum 
tunneling of the variable 9 between different topological 
sectors, e.g. between the states 9 = and 9 = 2ir, are 
defined as solutions of the equation SS/S9 = 0. For the 
action (Q), (f7|) this equation is approximately satisfied 
for a wide class of sufficiently smooth functions obey- 
ing the boundary conditions 9(0) = and 9(/3) = 2tt. 
Defining the typical instanton frequency Q, one finds 
§,0 that for Q, < E c and T <C E c the total ac- 
tion S[9(t)] does not depend on Ec and is equal to 
Sinter)] = Airar, while for f2 > Ec the instanton con- 
tribution gets suppressed by the kinetic energy term Sq. 
For Si n t[9(T)] ^> 1 both the tunneling amplitude between 
the states 9 = 2nm and 9 = 2ir(m ± 1) and the flux- 
dependent part of the free energy F are exponentially 
reduced as cx exp(— S[9(t)]). Hence, in the limit of small 
4> x <C 1 we obtain / = eE c 4> x /TT, where 



E c /E c = A(T) exp(-47rar). 



(13) 



This result is justified provided Altar 3> 1 and the ratio 
E* c jEc is sufficiently small. It demonstrates that also 
at low T <C Ec PC is strongly suppressed by Coulomb 
interaction provided the ring perimeter 2ttR exceeds the 
dephasing length L v ([l2|). Eq. ( |l3| ) suggests that at low 
T interaction-induced decoherence in our model is con- 
trolled by the parameter ar ~ a Y] 1 na n rather than 
by a or odnr. Hence, for r>l the decoherence effect of 
Coulomb interactions is much stronger than it could be 
expected both from perturbative and RG approaches. 

While the above instanton analysis is applicable for 
Anar 3> 1, it can hardly provide a quantitative descrip- 
tion of PC for moderate ar. Furthermore, even for large 
rings rigorous analytic evaluation of the pre-exponent 
A(T) in eq. ( [i"3| ) is difficult because the action ([j]), (Q) is 
strongly non-Gaussian and the standard techniques are 
hard to employ. On a qualitative level one expects that 
many non-Gaussian quasi-zero modes should yield a large 
pre-exponent A(T) increasing upon decreasing T. 

Quantum Monte Carlo analysis. In order to quanti- 
tatively describe reduction of PC at arbitrary values of 
ar we carried out extensive Monte Carlo (MC) simula- 
tions and directly computed E c as a function of a and 
r. Our method is described in details in ref. jO]. In 
order to discretize the action defined by eqs. (||), (||) 
and (0) we introduce the Trotter number TV = (3Ec- We 
then evaluate path integrals in eq. (Q) for different m 
and determine E C (T) = 27r 2 T(m 2 )0 x= o. By increasing 
N we decrease effective temperature in our simulations. 
The value E C (T) grows with decreasing T and finally 
reaches a plateau which defines the zero temperature 
value E' c (a,r), see the inset in Fig. 1. 
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Fig. 1 The value Eq/Ec as a function of r for a — 0.019 
and different T. From bottom to top: T/E c = 1, 0.5, 0.14, 
0.06, 0.03 and 0. The inset: E* c /E c as a function of T. From 
top to bottom: r = 10, 30, 60, and 120. 
At relatively high temperatures T > arEc our MC data 
confirm the semiclassical result ( |ll| ) . Also at lower T we 
observe dramatic suppression of E c with increasing r, see 
Fig. 1. At sufficiently large r this suppression turns ex- 
ponential and at Ec / (47rar) < T < Ec is well described 
by eq. @ with A(T) w exp(c£;c/T) and c ~ 1. Thus, 
eq. (|l^) holds even at temperatures well below the level 
spacing (~ Ec) of the noninteracting problem. 
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Fig. 2 Eq/Ec as a function of r at T = (bottom to top: 
a = 0.21, 0.07, 0.019, and 0.005). The inset: the same data 
for different a collapse onto one curve if plotted versus ar. 
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The data corresponding to the exact limit T = are 
displayed in Fig. 2 for different a. At relatively small 
r we observe very weak suppression of Eq in agreement 
with our perturbative result (|^). However, at ar > 1 the 
ratio E* c jEc rapidly decreases as Eq/Ec oc r~ 7 with 
7 rs 1.8. In order to judge whether at T = such power 
law behavior holds at all r > 1 /a or it is merely an inter- 
mediate regime with further crossover to the exponential 
dependence ( |l3| ) (cf. Fig. 1) it is necessary to extend 
our simulations at large N to values Eq/Eq below 0.01. 
Such details are, however, of no significance for our con- 
clusions. Much more importantly, our MC data clearly 
demonstrate drastic suppression of PC by Coulomb in- 
teractions even exactly at T = in rings with 2nR > L v 
as well as the failure of perturbation theory for such val- 



ues of R. Hence, L v (|12|) remains the relevant scale down 
to T — 0, see the inset in Fig. 2. 

Real time calculation. In order to illustrate a rela- 
tion between the above calculation and real time qua- 
siclassical analysis let us evaluate the return probabil- 
ity for a pair of classical paths which encircle the ring 
and return to the initial point after time t. As in Q 
we distinguish two probabilities W\ and W 2 correspond- 
ing respectively to two identical and two time-reversed 
paths. The quantity W\ is the standard return prob- 
ability and W 2 describes quantum interference between 
time-reversed paths. The latter quantity is sensitive to 
quantum coherence and it vanishes in the classical limit. 

Without interactions the time reversal symmetry 
yields W\ = W% — W^-°\ The effect of interactions 
is accounted for by the influence functional |l3| T = 
exp(iSfi ~ Si), and the probabilities are now equal to 
W(°'J-. For our model the actions Sr and Si are eval- 
uated by integrating out the fluctuating electromagnetic 
field produced by the electron gas. Similarly to eqs. (Q), 
(||) both Sr and Si are expressed in terms of the inverse 
dielectric function l/e(u>, k) of the electron environment. 

For two identical paths one finds Sr = Si = 0, i.e. the 
probability W\ — is not affected by interactions 

at all. For a pair of time-reversed paths one also has 
Sr, = 0, but the action Si is now positive Si > and, 
hence, W% — exp(— In order to evaluate Si it 
is convenient to "unfold" the time-reversed path into two 
straight lines f i(ii) = vt\ and ^2(^1) = v(t — ti), where 
< t\ < t. Inserting these paths into Si S we get 
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where! = 2 cos(kvt-) — 2 cos[kv(t— 1 + )} &ndt± = ti±t2- 
In the long time limit one finds 
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Defining the path length L — vt, at T — we obtain 
e 2 L 



Si = 



I6ir 2 al 2 



L 
la 



(16) 



Hence, similarly to PC, at T — > the probability W2 is 
suppressed by interactions beyond the length ( |l2[ ) and 
the ratio W2/W1 oc exp(— L/L v ) vanishes at large L. 
The latter result proves that reduction of PC is due to 
zero temperature decoherence and is not just a renormal- 
ization effect. The same conclusion follows from studying 
the AB conductance which we do not present here. 

Finally, let us define the dephasing time t v . For bal- 
listic particle motion considered here we set t v — L v /v. 
In particular, for v ~ vp we get t v ~ T e (kpl) 2 , where 
r e = 1/vf is the elastic scattering time for electrons. The 
latter expression for t v is identical to the electron dephas- 
ing time obtained in ref. for 3d disordered conductors. 
This observation is consistent with the idea |l4| that the 
results of weak localization and persistent current exper- 
iments may be closely related. In particular, we expect 
both to yield the same dephasing length L v at low T. 

We are grateful to F. Guinea, M. Paalanen and G. 
Schon for interesting discussions. 
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